Spin wave dynamics and the determination of intrinsic damping in locally-excited 

Permalloy thin films 
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Time-resolved scanning Kerr effect microscopy has been used to study magnetization dynam- 
ics in Permalloy thin films excited by transient magnetic pulses generated by a micrometer-scale 
transmission line structure. The results are consistent with magnetostatic spin wave theory and 
are supported by micromagnetic simulations. Magnetostatic volume and surface spin waves are 
measured for the same specimen using different bias field orientations and can be accurately calcu- 
lated by fe-space integrations over all excited plane wave components. A single damping constant 
of Gilbert form is sufficient to describe both scenarios. The nonuniform pulsed field plays a key 
role in the spin wave dynamics, with its Fourier transform serving as a weighting function for the 
participating modes. The intrinsic Gilbert damping parameter a is most conveniently measured 
when the spin waves are effectively stationary. 

PACS numbers: 75.30.Ds, 75.50.Bb, 75.70.Ak, 76.50.+g 



Interest has been growing for many years in the time- 
domain investigation of magnetization dynamics in re- 
sponse to a short magnetic pulse. An impulse excitation 
is broadband, but the dynamics of the magnetic system 
are also very sensitive to the parameters of the excitation 
pulse. Time-domain pulse shaping can eliminate ring- 
ing for a coherently-switched ferromagnetic element, to 
greatly enhance the performance for applications [H, Q . 
In addition, the spatial profile of the pulsed field is critical 
in dictating the magnetic dynamics. Magnetostatic spin 
waves generated by such nonuniform transient field have 
been observed in a number of time-resolved optical [3| 
and inductive @, @] experiments, and also in frequency- 
domain studies [6]. The resulting position-dependent 
temporal response creates additional challenges for char- 
acterizing the dynamics and for determining the intrinsic 
magnetic damping. The focus of the present work is to 
address these difficulties within a simple physical frame- 
work. 

Under the condition of linear behavior of the spin wave 
dynamics in the low amplitude regime 0,0], the magnetic 
response is the linear superposition of all plane wave com- 
ponents that can be excited by the pulse. Assuming the 
spin waves propagate only along the x direction (the co- 
ordinate system is defined in Figfjja)), the out-of-plane 
component of magnetization can be described by: 



M z (x,t) = e 



P{k,uj(k))sin[kx-uj(k)t+(j)]dk (1) 



where r _1 is the inverse decay time, ui(k) is the disper- 
sion relation, and P(k,uj(k)) is the spectral density of 
the pulse field determined from its spatial and temporal 
profiles and acts as a relative weighting factor for the dif- 
ferent spin wave components in fc-space. The influences 
of the pulse field parameters, the intrinsic damping, and 
the dispersion relation of the spin waves are contained 
explicitly. k c is a cut-off wave number for numerical inte- 
gration (fc c = 5 /im" 1 is sufficient for the magnetostatic 



regime with the stripline dimensions used here). The 
initial phase angle <j> is taken to be independent of fre- 
quency on account of the pulse excitation. The longer 
trailing edge of the pulse causes a non-oscillatory quasi- 
static background that is not included in Eq.(TT]), but in 
actual calculations we take the pulse shape into account 
by fitting the high-frequency data Q. 

The experimental geometry is shown by the inset 
schematics in FigfT] A Ni§oFe2o film (Permalloy, or 
"Py"), with thickness d = 10±1 nm, is deposited on a 150 
fim thick glass substrate using e-beam evaporation. The 
film is then clamped on a copper thin film stripline, with 
a small amount of high- vacuum dielectric grease applied 
for a strong surface-tension bond and to ensure electri- 
cal isolation between Py and Cu. Two coplanar wires in 
the stripline structure transmit a current pulse (rise time 
< 20 ps at the sample) from a semi-insulating GaAs pho- 
toconductive switch (carrier lifetime ~ 300 ps) and gen- 
erate a nonuniform magnetic pulse h(x,t). The width 
and separation of the wires are both 3 fim and are much 
smaller than the length (~ 400 //m), and the system 
can be treated as quasi-one-dimensional. The pulse field 
amplitude (and corresponding initial torque) decreases 
quickly away from the wires, falling to less than 10% of 
the peak value beyond \x\ = 8 /jm for the experimental 
geometry. We define the region of the Py film enclosed 
by these boundaries as the "source" area. An in-plane 
bias field Hq saturates the magnetization of the Py film 
such that Mq is parallel to Hq. This layout fixes the 
directions of the wave vectors k to be parallel to rc-axis, 
and enables us to detect different spin wave modes by 
changing the direction of H a (M ). The focus here is on 
the special cases of k \\ M Q and k _L Mo, although other 
angles can be similarly addressed @] . Changes of M z are 
measured by means of time-resolved scanning Kerr ef- 
fect microscopy Q. This technique offers ~ 500 nm spa- 
tial resolution determined by the spot size of the focused 
probe beam (much smaller than typical spin wave length 
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encountered in this work), and introduces a versatility 
relative to the pulsed inductive method, whose spatial 
resolution is limited by fixed-position, micrometer-size 
probe devices, and suffers a loss of signal when the mag- 
netization is perpendicular to the wires. The optical ap- 
proach allows a direct determination of a variety of spin 
wave dispersion laws. 

Quasi-ID micromagnetic modeling was also carried out 
in order to benchmark the /c-space calculation. The mag- 
netic film was discretized along the x-direction, such 
that the "finite elements" were 10 nm in both x- and 
z-directions, while infinitely long in y-direction. The 
spin motion of each element obeys the Landau-Lifshitz- 
Gilbert (LLG) equation [loj]: 

dM 



tt a dM 

- 7o M x H c ff + — M x — — 

M, at 



(2) 



where 70 = 17.61 MHz/Oe is the gyromagnetic ratio, 
M s = 760 emu/cm 3 is the saturation magnetization of 
the Permalloy film, and a is the Gilbert damping pa- 
rameter. H c ff is the effective field mainly contributed 
by the external field and the magnetostatic field. The 
exchange interaction is found to be insignificant in the 
magnetostatic regime [3, dH . The implementation of the 
simulation follows standard procedures 13, 12 1. 

Figfjja) shows the typical response of M z , in the 
k || Mo geometry, with Hq x = 200 Oe. The solid curve 
is a measurement at x = 0, and other positions show al- 
most the same profile but with amplitude decrease with 
increasing x. The damping is spatially uniform and not 
influenced by the propagation of spin waves [l3L 14 1 , con- 
sistent with expectation for k || M$. This geometry 
exhibits magnetostatic backward volume waves with a 
dispersion law (l5j : 



ujj +LU H LV M (l-e~ kd )/kd 



(3) 



where ujh — JqHq, u>m = ^loMs- For the present ex- 
periments, the group velocity, v g = is on the order 
of 0.1 /xm/ns and the spin waves are effectively station- 
ary over the time scale of the measurement. The cal- 
culated waveform based upon Eq.([T]) is plotted by the 
crosses in FigQTa), ancl agrees well with the measured 
data. To determine the weighting function P(k, ui(k\), a 
standard linearization analysis is applied on Eq.([2]) [161 ] 
to give M z = ffi^ h(u)h z (k), where h(uj) and h z (k) 
are the Fourier transformation of h(t) and h z (x), respec- 
tively. The w-dependent quantities are found to be in- 
significant in the calculations (except for a negative sign 
if lj > ujh), so the weighting function can be approxi- 
mated by P(k) = \h z (k)\. The Biot-Savart law was used 
to calculate the in-plane (h x ) and out-of-plane (h z ) com- 
ponents of the excitation field. The spatial distributions 
of h x and h z depend on the distance A between the plane 
of the Py film and the plane of the stripline. A cannot be 
precisely measured here and is used as a fitting param- 
eter. In this nearly stationary case, magnetic damping 
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FIG. 1: Different damping behavior of local magnetization 
under nonuniform excitation. The coordinate system is de- 
fined in (a) and is the same throughout the work. The paral- 
lel rectangular bars represent the stripline structure and the 
black dots represent the probe points (not to scale). Rela- 
tion between the wave vector and the bias field is shown in 
each panel. The solid curves are measured M z (t) traces, the 
crosses are calculated results based on Eq.([T]), and the open 
circles are simulated results using the quasi-lD LLG model 



of the system is unambiguously determined by the ex- 
ponential decay time r in Eq.([T]), and can be measured 
directly from the logarithm of the decreasing amplitude 
of the experimental waveform; the result for Fig[TJa) is 
t = 1.40 ± 0.01 ns. On the other hand, the Gilbert 
damping parameter a can be independently fitted by mi- 
cromagnetic simulations based solely on LLG equation; 
the result for this case is a = 0.0081±0.0003, and the sim- 
ulated waveform is plotted by open circles in FigQJa), in 
excellent correspondence with the measurement and the 
/c-space calculation. The two damping parameters are 
related by r = [a"/o(2iTM s + Ho)] -1 [111, and our results 
obtained by independent fittings are consistent with this 
relation. 

Good agreement between the measurements, fc-space 
calculations, and quasi- ID simulations are again obtained 
in the k _L Mo geometry, as shown in Figfjjb-d) f° r the 
probe positioned at x = and Ho y ranging from 110 Oe 
to 300 Oe. In this geometry magnetostatic surface waves 
(MSSW) are the dominant modes, leading to qualita- 
tively different spatiotemporal dynamics in the Py film. 
For the /c-space calculations, the dispersion law of MSSW 

lo 2 = {uj h + ujm/2) 2 - {uj A4 /2) 2 exp(-2fcd) (4) 

is used to calculate M z (x,t) using Eq.|T|). Linear 
response theory for the surface modes yields M z = 
oj'^-u 2 ^ l ( UJ )[ UJ Hh z (k) + iuih x (k)], that is, the out-of-plane 
magnetization responds to both in-plane and out-of-plane 
pulse fields. In the present work u>h <C ui, the contri- 
bution of h z is small and the weighting function again 



can be approximated by P(k) = \h x (k)\ [2(| (the < 
dependence is neglected as in the k || Mo case). Tl 
M z (t) traces show significantly shortened decay tin 
(compare FigQJb) to FigOJa)), since the excited surfa 
modes possess fairly large group velocity to transfer tl 
noncquilibrium spin wave energy out of the probed pos 
tion. At x = 0, the decay time shown in FigQJb-d) do 
not change explicitly when Ho y decreases (which leads 
increasing group velocity), but after we average M z (x, 
over a; | < 8 fj,m, the decay time in the whole "souro 
area indeed decreases with larger group velocity, as e 
pected (results not shown). In other words, the dam 
ing behavior in MSSW configuration cannot be quantit 
tively described by single-point measurements, but has 
be analyzed using the global approaches (scanning prol 
experiment, micromagnetic simulation and /c-space calc 
lation). The damping behavior is naturally embedded 
Eq.([T]), because of dephasing of the component freque 
cies through the u>(k)t term. Eq.JT]) is a generalization 
the formula proposed in Ref.[4j], which gives an intuiti 
description of the phenomenon. The damping envelope 
is jointly determined by a Gaussian term accounting for 
spin wave dispersion, and the intrinsic exponential de- 
cay. Our general approach works for both k \\ Mq and 
k _L Mq geometries, and explains why the decay remains 
exponential when the spin wave propagation is negligible. 

In the k _L Mq geometry, individual wave packets 
propagating in the x direction can be observed when 
M z (t) is measured outside the "source" area (i.e., the 
probe point moves away from the stripline). Represen- 
tative results for the case of Hg y = 80 Oe are shown in 
FigI2a-d). Recording a two dimensional (position and 
time) map for the peak of the wave packet, its group 
velocity can be determined to be 4.8 ± 0.5 fim/ns. Per- 
forming such analysis for a range of bias fields yields the 
group velocity as a function of frequency, as shown in 
FiglUe). The measurements and numerical calculations 
agree reasonably well with the MSSW theory (v g deter- 
mined from Eq.Q). 

For the propagating wave packet discussed above, the 
M z (t) is asymmetric in time (the increase of the os- 
cillation amplitude appears "slower" than the follow- 
ing decline). This asymmetry is especially apparent 
in Fig[2^c,d), but again is well reproduced by the k- 
space calculation based on Eq.([T]) (and which cannot be 
achieved by the Gaussian- type formula in Ref.[4]). This is 
shown by the crosses in Figf2ja-d), and is also supported 
by the quasi-lD simulations (open circles). A single scale 
factor fits the measured amplitude at all positions in the 
considered range, indicating that the Gilbert mechanism 
(a = 0.0081 is still used here [2lj]) also accounts for the 
spin wave attenuation during propagation 0, [l3| . Mea- 
surements of this attenuation are potentially an effective 
way to determine a when k _L Mq. 

The temporal profile of M z (t) ultimately stems from 
the spatial profile of the excitation field, which directly 
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FIG. 2: Propagation of a spin wave packet in the k _L Mo 
geometry, (a) - (d): M z {i) traces at x = 10, 12.5, 15 and 
17.5 /xm, respectively, when Ho y — 80 Oe. The solid curves 
represent the measured data, which are extracted from a sin- 
gle spatiotemporal scan, the crosses are the calculated results 
based on Eq.fl}, an d the open circles are the simulated re- 
sults using the quasi- ID LLG model. The simulated and cal- 
culated waveforms are scaled to the same amplitude as the 
experimental data, (e): Group velocity of MSSW modes as a 
function of oscillation frequency, with Ho y ranging from 30 to 
300 Oe. The black dots are measured data, and the crosses 
are results from quasi-lD simulations. The grey shadow is 
the theoretical curve calculated from MSSW dispersion law 
(Eq.Q), with the width reflecting its lower and upper lim- 
its due to the uncertainty in the cutoff wavelength and film 
thickness (k — ~ 5 /im" 1 and d = 9 ~ 11 nm). 

determines P(k) in Eq.fl]). Fig[3] illustrates the effect of 
different distributions P(k), calculated for several values 
of the film-stripline gap A. As A increases, the spatial 
variation of the pulse field becomes smoother and P(k) 
acquires relatively higher spectral density at smaller k 
(FigUHa)). This yields relatively larger oscillation ampli- 
tude at early times before the wave packet peaks, and can 
be understood here as a consequence of more spin wave 
components with higher phase velocity v p = FigsGJb- 
d) show the bracketing of the "best-fit", A = 1.6 fim, to 
the experimental data as presented in Fig(2)(d). 

As the spin wave packet propagates farther away from 
the source, its shape is gradually broadened and the peak 
time eventually exceeds the maximum optical delay of 
the apparatus (5 ns). In this situation, the oscillations 
are dominated by small-fc components. The experiments 
yield a unique position at each bias field (x ~ 30 /mi for 
the case of Hq v = 80 Oe) where the incoming energy from 
the propagating mode effectively balances the intrinsic 
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FIG. 3: Influence of the spatial distribution of the pulse field, 
(a), Calculated distributions of normalized P(k) for A = 0.5 
/im (dotted curve), A = 1.6 /im (solid curve), and A = 2.5 
/im (dashed curve), (b) - (d), M z (t) traces at x — 17.5 /im 
calculated with Eq.JTJ, using A = 0.5 /im, A = 1.6 /im, and 
A = 2.5 /im, respectively. 



dissipation at that position, such that a nearly time- 
independent oscillation amplitude is observed through- 
out the measurement window, as shown in FigJUa). For 
larger x, the power balance is broken and the intrin- 
sic decay of the long-wavelength oscillations dominates 
(Fig0[b)). 

In summary, we have studied the spatiotemporal dy- 
namics of a ferromagnetic film in response to a short 
magnetic pulse localized in one spatial dimension. The 
response can be interpreted as a superposition of plane 
waves modulated by the spectral densities of the spatially 
nonuniform, transient excitation field, and decaying ac- 
cording to intrinsic Gilbert damping. Magnetostatic vol- 
ume modes and surface modes can be self-consistently de- 
scribed within this interpretation, and the experimental 
and analytical results show good agreement also with mi- 
cromagnetic simulations. The Gilbert damping param- 
eter could be directly measured from broadband FMR 
waveforms when the spin waves are effectively stationary. 
The fc-space calculations also require negligible computer 
resources in comparison to micromagnetic simulations, 
and offer an opportunity to invert the problem and de- 
sign magnetic waveforms for applications. In addition, 
the analysis can be extended to a second dimension to 
account for thicker or multilayered structures, or for 2D 
spin wave propagation. 

This work was supported by NSERC, iCORE, CIAR 
and CRC. The samples were fabricated in the University 
of Alberta Nanofab. 




2 3 4 
f(ns) 



FIG. 4: Spin wave oscillations measured at (a), x ~ 30 /mi 
and (b), x ~ 35 /im. The experimental conditions are the 
same as those in Fig[3a-d). 
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